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PREFACE 


This report presents the results of studies conducted during 
the period June 20, 1968 - July 19, 1969, under NASA research contract 
NAS 8-21432, "Lunar Surface Engineering Properties Experiment Defini- 
tion." This study was sponsored by the Advanced Lunar Missions 
Directorate, NASA Headquarters, and was under the technical cogni- 
zance of Dr. N. C. Costes, Space Science Laboratory, George C. Marshall 
Space Flight Center. 

The report reflects the combined effort of four faculty investi- 
gators, a research engineer, a project manager, and six graduate 
research assistants, representing several engineering and scientific 
disciplines pertinent to the study of lunar surface material properties. 
James K. Mitchell, Professor of Civil Engineering, served as Principal 
Investigator and was responsible for those phases of the work con- 
cerned with problems relating to the engineering properties of lunar 
soils and lunar soil mechanics. Co- investigators were William N. 
Houston, Assistant Professor of Civil Engineering, who was concerned 
with problems relating to the engineering properties of lunar soils; 
Richard E. Goodman, Associate Professor of Geological Engineering, 
who was concerned with the engineering geology and rock mechanics 
aspects of the lunar surface; and Paul A. Witherspoon, Professor of 
Geological Engineering, who conducted studies related to thermal 
and permeability measurements on the lunar surface. Dr. Karel Drozd, 
Assistant Research Engineer, performed laboratory tests and analyses 
pertinent to the development of a borehole probe for determination 
of the in-situ characteristics of lunar soils and rocks. John 
Hovland, David Katz, Laith I. Namiq, James B. Thompson, Tran K. Van, 
and Ted S. Vinson served as Graduate Research Assistants and carried 
out many of the studies leading to the results presented in this 
report. Francois Heuze, Assistant Specialist, served as project 
manager and contributed to studies concerned with lunar rock mechanics. 



Ultimate objectives of this project are: 

1) Assessment of lunar soil and rock property data using 
information obtained from Lunar Orb iter and Surveyor 
missions . 

2) Recommendation of both simple and sophisticated in-situ 
testing techniques that would allow determination of 
engineering properties of lunar surface materials. 

3) Determination of the influence of variations in lunar 
surface conditions on the performance parameters of a 
lunar roving vehicle. 

4) Development of simple means for determining the fluid 
and thermal conductivity properties of lunar surface 
materials . 

5) Development of stabilization techniques for use in loose, 
unconsolidated lunar surface materials to improve the 
performance of such materials in lunar engineering 
application. 

The scope of specific studies conducted in satisfaction of these 
objectives is indicated by the following list of contents from the 
Detailed Final, Report which is presented in four volumes. The names 
of the investigators associated with each phase of the work are 
indicated. 


VOLUME I i 

MECHANICS AND STABILIZATION OF LUNAR SOILS 

1. Lunar Soil Simulation 

(W. N. Houston, L. I. Namiq, and J. K. Mitchell) 

2. Lunar Surface Traf ficability Studies 
(J. B. Thompson and J. K. Mitchell) 


3 . 


Foamed Plastic Chemical Systems for Lunar Soil Stabilization 
Applications 

(T. S. Vinson and J. K. Mitchell) 



VOLUME II 

LUNAR SOIL PROPERTIES FROM PHOTOGRAPHIC RECORDS 


v 


1. Soil Property Evaluations From Boulder Tracks on the Lunar 
Surface 

(H. J. Ho viand and J. K. Mitchell) 

2. Deduction of Lunar Surface Material Strength Parameters from 
Lunar Slope Failures Caused by Impact Events - Feasibility 
Study 

(T. S. Vinson and J. K. Mitchell) 

VOLUME III 
BOREHOLE PROBES 

1. The Mechanism of Failure in a Borehole in Soils or Rocks 
by Jack Plate Loading 

(T. K. Van and R. E. Goodman) 

2. Experimental Work Related to Borehole Jack Probe and Testing 
(K. Drozd and R. E. Goodman) 

3. Borehole Jack Tests in Jointed Rock - Joint Perturbation and 
No Tension Finite Element Solution 

(F. E. Heuze, R. E. Goodman, and A. Bomstein) 

VpLUME IV 

FLUID CONDUCTIVITY OF LUNAR SURFACE MATERIALS 

1„ Studies on Fluid Conductivity of Lunar Surface Materials 
(D. F. Katz, P. A. Witherspoon, and D. R. Willis) 



vi 


VOLUME IV OF FOUR CONTENTS 

Page 

Chapter 1. Studies on Conductivity of Lunar Surface 

Materials 1-1 

D. F. Katz , p. A. Witherspoon , and D. R. Willis 

I. Objectives 1-1 

II. Introduction 1-1 

A . Conceptual Description of Probe 1-1 

B. Flow in Porous Media in the Lunar 

Environment 1-3 

III. The Local Similarity Method 1-5 

A. Basic Theory 1-5 

B. Method for Determining Permeability and 

Area Fraction 1-8 

C. Discussion of Assumptions 1-10 

IV. Comparison with Existing Theory 1-15 

V. Proposed Experimentation 1-20 

VI. Conclusions 1-22 

References 1-24 

Appendix A — Average Quantities for a 

Symmetric Flow 1-25 

Appendix B — Knud sen Number 1-27 

Symbols 


1-28 



CHAPTER 1 


1-1 


STUDIES ON CONDUCTIVITY OF LUNAR SURFACE MATERIALS 
(D. F. Katz, P. A, Witherspoon, and D. R, Willis) 

I. OBJECTIVES 

The overall objective of this investigation is to develop a 
means of measuring the permeability of lunar soils and rocks in situ. 

It has been proposed to design and test a surface probe that can 
measure permeabilities with reasonable accuracy. Because of the lack 
of any atmosphere on the moon, it will be necessary to utilize gas 
in operating the probe. In the current year's work, a theory of gas 
flow in porous media appropriate to the surface probe has been 
developed. The theory is applicable to probe operation in the lunar 
environment as well as on earth. In addition, the design of experiments 
to check the theory has begun. 

II. INTRODUCTION 

A. Conceptual Description of Probe 

A schematic drawing of the probe is given in Figure 1-1. 

The system basically consists of a holding chamber containing 
pressurized fluid, a disc-shaped source, and pressure measuring devices 
imbedded along with the source in a circular, impermeable skirt. 

The holding chamber contains the charge of gas for an individual 
measurement, and is connected by valves to the disc source and to a 
larger gas storage tank. A single switch releases the charge of 
gas, and activates a timer linked to the surface pressure tap. Constant 
source pressure is achieved by a servo connection between the release 
valve and an exit tap. The connecting valve between the holding 
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chamber and the storage tank insures that only a small portion of 
the gas supply is consumed per measurement. This also allows for 
variation in source pressure as required by varying permeability. One 
can then compute permeability from a record of pressures and flow rates. 

In the design of the probe, due consideration must be given 
to size, weight, durability requirements, etc. Conceptually, it 
seems feasible to utilize a waste gas, such as C0 2 « In view of current 
practice on earth, a miniaturization of the probe is suggested. 

However, any miniaturization is limited by the requirement that the 
source diameter be large compared to the pore dimensions. 

B. Flow in Porous Media in the Lunar Environment 

Due to the high vacuum in the lunar environment, there exists 
a basic problem of determining the dominant flow regime during probe 
operation. As the average pore size and/or fluid pressure decrease, 
the fundamental nature of gas flow changes. The initial departure 
from viscous flow is the relaxation of the no slip boundary condition 
on the interstitial surfaces. The resulting, augmented flow, while 
still viscous in nature , exhibits a greater temperature dependence 
than previously. As the effective degree of rarefaction increases 
further, the continuum nature of the fluid breaks down, and the flow 
must be modeled from a molecular point of view. In general, account 
must be taken of both intermolecular collisions, and collisions 
between molecules and the interstitial boundaries. However, when a 
high degree of rarefaction is achieved, the former become negligible 
due to the relative scarcity of fluid molecules. The fluid parameter 



indicative of the particular flow regime prevalent is the Knudsen number 

defined here as Kn_ 5 where X is the mean free path of the molecules 

d d 

of the fluid, and d is an average pore cross sectional dimension. For 
viscous flow, Kn^ « 1, whereas for the so-called Knudsen flow, Kn ^ » 1 

The geometrical complexity of the pore structure of rocks 
precludes "exact" solutions of the fluid equations of motion, and one 
is forced to work in terms of averaged flow quantities; see Appendix A. 
This approach has proven successful in treating terrestrial, continuum 
flows. It is, therefore, not unreasonable to expect that transitional 
flow in the lunar environment may also be amenable to such an approach. 
Henceforth, all flow quantities referred to are averages. It is also 
necessary to assume that the porous medium of interest is homogeneous 
and isotropic. 

It seems likely that in the immediate neighborhood of the 
fluid source, the flow will be viscous in nature, and, thus, basically 
predictable by current techniques. However, the nature of the 
evolution of Knudsen flow, as the distance from the source increases, 
is extremely difficult to foresee at this time. Detailed theories 
for transition flow, namely solutions of the Boltzmann equation, for 
even the simplest geometries are quite scarce. As a result, this 
report introduces the concept of local similarity in treating the gas 
flow. This approach makes maximum use of physical intuition about the 
flow field, and requires experimental calibration of the probe on a 
sample of known material. It does not, however, require actual solution 
of the fluid equations of motion. Such an approach to the problem 
seems far more reasonable, as a first effort, than some approximate 


over-sinplified solution. 
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III* THE LOCAL SIMILARITY METHOD 
A. Basic Theory 

In the analysis of terrestrial fluid flow in porous media , 
effective use if made of a simplified, integrated form of the equation 
of motion in which the velocity, or flow rate, is taken to be 
proportional to the pressure gradient. For continuum flow this is 
known as Darcy* s Law, and can be written 

v = C(p)Vp (1-1) 

where v is the velocity vector and p is pressure. The coefficient C 
contains parameters that characterize the particular problem, basically 
viscosity and appropriate pore dimensions. For incompressible, continuum 
flow, C is independent of pressure, and is normally written 

C = - (1-2) 

U V 

where y is viscosity, and k is termed the permeability. Thus, 
permeability depends upon the size and shape of the interstices, and 
has dimensions of length squared. For both slip flow and Knudsen flow, 

C is pressure dependent (Carmen, 1956). It is, therefore, suggested 
that Equation (1-1) may be characteristic of flow in porous media in 
general. This assumption is made, and thus C(p) can be thought of 
as a "master 11 diffusion coefficient encompassing all flow regimes. 

C (p) can be rendered dimensionless by defining 
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C(p) = C (p) (1-3) 

where C is a dimensionless function of pressure. It is assumed that 
C is a function of local pressure (or Knudsen number) * only. L 2 has 
dimensions of area, and can be thought of as a generalized permeability, 
applicable to all flow regimes. These are the basic local similarity 
assumptions. They embody all the necessary physics of the flow field. 
The validity of these assumptions will be discussed below. 

Consider now a steady, isothermal flow field dependent upon only 
one space variable, denoted by r, i.e. a one-dimensional, cylindrically 
symmetric, or spherically symmetric problem. Henceforth such problems 
will be referred to as "symmetric." Equation (1-1) becomes 

v = C(p)|£ (1-4) 

The continuity equation can be written 

apvr-* = Q (1-5) 

where a is the area fraction, of. Appendix A j p is the average density 

✓s. 

and Q takes on different values (see Table 1-1) , depending upon 
the nature of the symmetric problem. 


^Pressure and Knudsen number are inversely proportional. See Appendix B . 
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TABLE 1-1 



Type of Symrnetric Problem 

Symbol 

One-Dimensional 

Cylindrical 

Spherical 

j 

0 

1 

i 

2 

Q 

Q/A r 

— ! 

Q/2^c 

Q/4ir 


Here Q refers to the dimensional mass flow rate, A , is the total 
x r 

cross-sectional area of the rock in one-dimensional flow, and % is 

c 

the cylinder length for cylindrically symmetric flow. Combining 
Equations (1?3) , (1-4)^ and (1-5). 


2 = 


a 


hi 

u 


pr j C(p)|2- 


(1-6) 


Assume that the gas is perfect, so that p = pRT, where R is the gas 
constant and T is the temperature. Then 


Q = 


a 


j/ 

yRT 


r j pC(p) 


d£ 

dr 


(1-7) 


Since C(p) is dimensionless, it can be expressed as a function of a 
dimensionless pressure £ s 


C(p) 


_ 1 

" F(0 


(1-8) 












The particular form of C is motivated by the fact that C can be 
considered as being only a function of an appropriate Knudsen number 
(see Appendix B ) . Thus 
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= LV _ 

2JJRT P “ Kn 

L 


(1-9) 


— i 8RT 

where V = / the mean thermal speed, and Kn is the Knudsen number 

L 

based on the "length" L, Kn = ~ . Substituting in Equation (1-7) and 

L L 


rearranging , 


aL 2 ^i— r j p §£ = F(C) (1-10) 

QJJRT ax 

This is a fundamental similarity relation. F(C) is a "universal" 
function, in that data from all symmetric problems, plotted according 
to Equation (1-10) fall on the same curve. It should be noted that F (£) 
is not the only possible universal curve. Multiplication by any real 
function of £ yields an equally universal curve. 


B. Method for Determining Permeability and Area Fraction 

By developing dimensionless relationships in this manner, 
Equation (1-10) can be used to determine both permeability L 2 and the area 
fraction a. For example, if the flow field for the lunar probe can 

A 

be approximated as spherically symmetric, then j = 2 and Q = Q/4? T. 

Equation (1-10) becomes 


air 


4TT 


^QyRT 


2 dp 
r d — 
p dr 


= F(S) 


(1-H) 
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Since F(£) presumably holds for any symmetric flow in porous media , 
this function can be uniquely determined by appropriate experimentation 
in the laboratory. 

To apply Equation (1-11) to an unknown rock sample, one must fix 
the flow rate Q and measure temperature T, from which y = y(T) is 
easily determined. If pressure taps are appropriately spaced in the 

d.D 

skirt of the probe (cf. Figure 1-1) then r 2 p-j£- can be measured at two or 
more different values of r. This enables one to determine two different 
values of the bracketed expression in Equation (1-11) which are sufficient 
to determine both L and a. 

A simple procedure illustrating the method of calculation 
can be outlined as follows. From Equations (1-9) through (1-11) 


In £ = In p + In 


fv 

L 2y: 


2yRT 


(l-12a) 


In F(C) 


In 

r . ^ 

2 dp 
r p 

+ In fotL 2 — — ; 


drj 

1 QURTJ 


- > 

, i r>tr\ d ln r2 P 
d In F(S) _ 1 _ dr 

d In ? d In p 


(l-12b) 


(l-12c) 


If follows that curves of F(C) vs ? and r 2 p ^vs p, plotted on 
identical log-log scales, differ only in the positions of their 
respective origins of coordinates. Thus, the two measured values of 
r 2 p ^ are plotted versus p on log-log paper, and F(£) vs £ is plotted 
on identical paper. The two plots are then placed on top of each 
other, and maneuvered, keeping respective axes mutually parallel, until 
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the two experimental points lie on the known curve . The coordinates 
of the origin of the experimental plot relative to the universal 
one then clearly yield both a and L, See Figure 1-2 « 

As drawn in Figure 1-2, F(£) will be a monatonic function of 
However, a problem of uniqueness in determining a and L will arise 
if there are regions where, on the log-log scale, F(C) is locally 
linear. If the experimental points correspond to such regions, then 
the position of the origin of the experimental plot is not uniquely 
determined. This difficulty can, principle, be circumvented, however, 
by redefinition of the universal curve such that it has adequate 
curvature throughout . 

C. Discussion of Assumptions 

The concept of local similarity contains the following basic 
assumptions : 

(1) flow is steady 

(2) flow is independent of the initial pressure in the 
porous medium 

(3) a single length characterizes the fluid conductivity 
of the porous medium 

(4) flow is isothermal. 

The first and second assumptions are not independent. This is shown 
qualitatively in Figure 1-3. Clearly the flow field in the entire 
rock is not steady, due to the presence of the pressure (and density) 
“wave front," i.e. the boundary between those portions of the rock 
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FIGURE 1-3 QUALITATIVE TIME DEPENDENCE OF FLOW FIELD FOR SURFACE PROBE 



that contain gas molecules from the source, and those that do not. 

The velocity of propagation of this front will depend in part upon 
the initial fluid pressure in the rock. Sufficiently close to the 
source, however, it is not unreasonable to presume that time 
variations in the flow field will be less significant, and subsequently, 
the influence of the initial pressure. Strictly speaking, then, the 
local similarity theory can only be applicable in this "quasi steady" 
region. It is assumed that this region is large enough to permit 
the necessary pressure measurements. 

The assumption of steady flow for operation of the surface probe 

can be examined quantitatively for two limiting cases: (a) slightly 

compressible, continuum flow, and (b) Knudsen flow. In both cases 

the governing equation becomes the diffusion equation (Carman, 

1956) . A solution to the equation for the precise geometry of 

the probe (without the assumption of spherically symmetric flow) has 

been given by Selim et al. (1963). Steady state flow is effectively 

reached for T = 10 3 , where T is a dimensionless time which can be written 
4Dt 

T = — y • Here t is time, "a" is the source radius, and D can be thought 

of as the diffusivity appropriate to the particular flow regime. For 

slightly compressible, continuum flow, D E , w here <f> is porosity and 

k is the continuum permeability mentioned earlier. For Knudsen flow, 

1 dV 

D = j , i.e,/ D is simply the Knudsen diffusivity divided by 
porosity. 

For continuum flow, consider a typical situation where (j) - lo” 1 , 
a - 10 1 cm, y - 10 1 poise, p ~ 10 6 dyne/cm 2 , and lo” 11 cm 2 < k < lo” 8 cm 2 

- 1 l 7 —2 

(Continuum flow isn't expected for k < 10 cm .) Then 10 sec < t < 10 
sec, so that steady state flow is feasible. 
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Now consider a typical Knudsen flow situation with the same 
values of <j> and "a" as above* Let T - 10 2 °K so that V ~ 10 4 cm/sec* 

„ g M g 2 o 

For most rocks, 10 cm < d < 10 cm* It follows that 10 < t < lCr 

sec. So, for the rarefied limit, as well as the continuum one, steady 
state flow appears feasible in probe operation. 

It is also instructive to utilize the Selim et al. solution 
in considering the assumption of spherically symmetric flow. Their 
results show that, at the surface of the porous medium, for r < 5.633a 
the steady flow field appears spherically symmetric. Thus, the surface 
pressure taps in the probe would have to be replaced in this region, 
a quite reasonable design requirement. The Selim et al. results also 
show that over practical time periods, there is essentially no flow 
at the rock surface for r < (10a) . Thus, for a = 0.1 cm, a reasonable 
radius for the impermeable skirt might be, say, 3 cm. This requirement 
is also compatible with the envisioned probe design. 

The problem of characteristic lengths is fundamental to flow 
in porous media, and, even for purely continuum flows, is as yet 
unresolved. In the local similarity analysis there can be formal 
dependence upon only a single characteristic dimension, L. Current 
theory on flow in porous media often introduces two characteristic 
dimensions. One is associated with effective pore diameter (denoted 
above by d) , and the other, with effective pore length. If the 
proposed experimentation reveals the necessity of introducing a second 
characteristic dimension, it will have to be done empirically. 
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IV. COMPARISON WITH EXISTING THEORY 

The most systematic approach to the problem of transitional flow in 
porous media has been made by Mason et al. (1961, 1962, 1963, 1964, 1967) 
in the form of the "dusty gas models". This model treats the overall 
transport of a single or multicomponent gas in a porous medium as a 
multiple diffusion problem. The pore structure is considered to be 
one of the diffusing species, but is constrained to remain stationary 
as a uniform distribution of "dust" particles. The starting point 
in the analysis is a "master" diffusion equation governing the effects 
of viscous transfer of momentum in the diffusion of a multicomponent 
gas mixture. This equation was derived using the advanced "thirteen 
moment" approximation of the kinetic theory of gases (Zhdanov, 1962). 

Only steady one«*dimensional flows are considered, and inertial effects in 
the continuum regime are neglected. Porosity and tortuosity effects are 
introduced empirically. For isothermal flow of a single component 
gas, the following equation is derived: 


Q = 


r 


aip + b 

•L+.£aE, 

< 

|1 + C 2 pJ 

4 


d£ 

dr 


(1-13) 


Here aj , b, ci, cz are constants, one of which must be determined 
experimentally. Since Q = praA, where A is the area of this fluid 
source. Equation (1-13) can be rewritten: 


v 


RT 

aA 


r 

f 


•V 

1 

aip + b 

1 + CiP 


P 


[1 + c 2 p_ 



w 


✓ 


d£ 

dr 


(1-14) 


Equation (1-14). is consistent with the local similarity expression. 
Equation (1-4) , for one dimensional flow with 
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C(p) 


RT 

OA 


1 

P 


(aip + b 


1 + c iP~ 

1 + CzP' 


(1-15) 


Equation (1-15) has the same quantitative pressure dependence as 
an equation obtained by Wakao et al* (1965) . For steady, one-dimensional 
isothermal flow that is undergoing transition in porous media* The result 
of Wakao et al . can be expressed as 


r r 2 p r 

°L 


c( P ) = -! 


ug 


D. 


J^R_ 


KA Kn_ p 
d ^o 
o 


2- + Kn , 
P,, d . 


1 + 


_I_ E_ 


Kn, p 
d *o 
o 



s 

D 


KA 


i + — — 

"il 


p o 

o 

j 


(1-16) 


The subscript o refers to reference conditions at the gas source. The 
first term in Equation (1-16) corresponds to simple Poiseuille flow. Wakao' s 
second term, to slip flow, and the third term, to Knudsen flow. Wakao* s 
equation contains no disposable constants, but its justification is 
largely heuristic, and is based on the single capillary model for flow 
in porous media. Similarly, the additional efforts, listed in Reference 8, 
are essentially heuristic. Unfortunately, experimental verification 
for all the above work has been confined to measurements on single 
capillaries or bundles of capillaries. 

The only experimentation on transitional flows in actual porous 
media has been the recent work of Huang and Ramsey (1968) . Their experiments 
consisted of steady, one -dimensional flows, in which large relative 
pressure differences between the ends of the rock samples were achieved. 
However, all pressures remained of the order of atmospheric, and rarefied 
effects were assumed to be present only by virtue of the small pore 
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sizes of the rock samples,, The data are in good agreement with the 
Wakao equation, cf. Equation (1-16) * However, it should be emphasized that 
the applicability of the data to the present problem is questionable 
due to the lack of in-vacuo conditions. 

The final work of interest is that of Sreekanth (1968) . He conducted 
experiments on transitional flows in single tubes in which large 
relative pressure drops were achieved, the tube exits being maintained 
at vacuum conditions. Three different tubes were used, with length-to- 
diameter ratios of 4,77, 9.66, and 15.03. The flows were steady and 
isothermal . Sreekanth * s data for pressure distribution agree extremely 
well with a continuum analysis, utilizing an approximate method of 
accounting for compressibility and assuming velocity slip (Ebert and 
Sparrow, 1965) . Such an analysis is consistent with the local similarity 
approach, and an effective C (p) can be deducted. 


_ A 4RT o 

C(p) vpr$T 


if T 

£o 

3 

'b i ReKn d X’ 
o 

+ 

f 1 

fo 

2 

Rex 



r -n 

1Tb 1 

7T 

p 

1 > 


4 

J 

P 


2 

L 4 


i p J 

8ReKn _ 
d 

v o> 

64ReKn^ 

d 

o 


(1-17) 


Here Re - ~ a Reynolds number based on mass flow rate; bi is a 
it yd J 

constant, of order one, contained in the slip boundary condition 

is the normal derivative of 


(v) = — biA 

wall 


9v 

. where 

fH 

3n 

V. > 

| wall 

3.n 

V ^ 


(wall 

velocity at the wall; and X is a weakly varying function of local 


Knudsen number, which is effectively a constant of order one. 
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Due to the fact that the Wakao et al« approach, cf* Equation (1-16) « 
which is intended for flow in porous media, is derived from a capillary 
model, it is instructive to apply their equation to the Sreekanth data. 

In so doing, it is necessary to select only those results with small 
values of Re and large values of Vp/p e because these are appropriate 
to transitional flow in porous media. Here Vp is the pressure drop 
along the tube, and p^ is the background pressure. If calculated 
volumetric flow rates V are compared with the experimental data, the 
accuracy improves with increasing length-to -diameter ratio. This is 
to be expected since the Wakao theory assumes fully developed flow, 
independent of capillary entrance effects. Typical results for the 
longest of the three tubes are shown in Table 1-2. As expected, the 


TABLE 1-2 


COMPARISON OF VOLUMETRIC FLOW RATES 
CALCULATED FROM WAKAO THEORY WITH SREEKANTH DATA 


P e 

* f 

V cm 

calculated min 

V / 

(Wakao et al .) 

v 

measured[rmn 

(Sreekanth) 

Percentage error 

% 

1.45 

78.76 

71.64 

9 

2.84 

639.52 

556.00 

13.1 

7.86 

1199.54 

990.50 

17.4 

10.25 

335.64 

293.00 

12.7 

20.05 

31.254 

29.48 

5.6 

26.53 

6.02 

5.95 

1.2 
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greatest accuracy is for the highest and lowest V (i*e» the lowest 

P e 

Reynolds number) * 

However , the calculated pressures along the tube length do not 

compare nearly as favorably with the data for large values of In 

p e 

some instances, these pressures become double -valued near the tube end. 
Typical results for the longest tube are shown in Figure 1-4. Clearly, 
then, the approaches of Sreekanth and Wakao et al . are not equivalent. 
Nonetheless, while there is some question as to which is more applicable 
to flow in porous media, neither contradicts the local similarity 
hypothesis made in this work. 



z 

J 


FIGURE 1-4 TYPICAL COMPARISON OF WAKAO AND SREEKANTH FORMULAS 
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V. PROPOSED EXPERIMENTATION 

In the final quarter of this year's work, plans were made for a 
series of experiments aimed at developing a probe prototype. To 
simulate the effects of the lunar environment, all experiments will be 
conducted on rock samples enclosed in a vacuum chamber. To our knowledge, 
this will be the first work of this kind that has been attempted. 

Facilities currently available in the Aeronautical Sciences Laboratories 

are easily able to maintain significant flow volumes at a vacuum of 10 6 

-8 

Torr , and with the aid of auxiliary equipment, a vacuum of 10 Torr 
can be reached. These pressure levels will enable us to examine the 
effects of rarefied gas flow without having to duplicate the much lower 

- 1 4 

vacuum (10 Torr ) of the lunar environment. 

The first experiment will consider one-dimensional flow. This is 
the simplest geometry, and will enable us to develop a familiarity with 
appropriate procedures for handling gas flow through porous media under 
high vacuum conditions. Two basic rock configurations are being considered. 
The first consists of a single, cylindrical rock specimen, with pressure 
taps embedded in the rock, cf. Figure l-5a. In the second configuration, 
the pressure taps are placed in gaps between cylindrical "wafers" of rock, 
cf . Figure l-5b, in order to insure that measurements of average pressure 
are obtained. 

By selecting rocks of appropriate permeability and controlling absolute 
pressures at the necessary levels, it should be possible to obtain conditions 
that vary from continuum viscous flow to flow with a high degree of 
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FIGURE 1-5 SCHEMATICS OF PROPOSED CONFIGURATIONS FOR 
EXPERIMENTS ON ONE-DIMENSIONAL FLOW 
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rarefaction « In so doing, the results obtained will provide a necessary 
test of the local similarity method* 

One of the fundamental considerations in the design of the lunar 
permeability probe is its range of application. Over what range of 
permeabilities can the lunar probe be expected to yield satisfactory 
results , given the operating conditions that prevail on the lunar surface? 
Can the probe be applied to consolidated and/or unconsolidated rocks? 

These problems can be investigated by varying the permeability of rock 
samples , and noting those conditions of temperature, flow rate, and 
absolute pressure that are practical for the envisioned probe. 

Another consideration is the type of gas to be used. Is there any 
advantage to be gained by using a single component gas rather than a 
gaseous mixture? Prom the experimental standpoint, it will be possible 
to use both approaches and thus provide data for another design factor 
of the probe. 

In subsequent experiments, the applicability of the envisioned probe 
to spherical flow conditions will be investigated. This work will provide 
design criteria for such features of probe construction as the necessary 
size of the probe skirt, and the optimum positioning of pressure taps. 

The results of this work can then be used to build and test a probe 
prototype. 

VI. CONCLUSIONS 

A theory of gas flow in porous media appropriate to the lunar surface 
permeability probe has been developed. This theory is based on the concept 
of local similarity and utilizes the assumptions of steady, isothermal, 
and symmetric flow assumptions that have proven possible for terrestrial 
flows. This concept of local similarity is more general than the few 



existing theories for transitional flow in porous media. It is entirely 
consistent with what is judged to be the best of such work. The design 
of experiments to check this theory and develop a probe has begun. 
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APPENDIX A. 


AVERAGE QUANTITIES FOR A SYMMETRIC FLOW 
In Equation (1-5) , the area fraction a is defined 


a 


A r. 

— P. 1 

A (rj) 


where 


A(rj) 


L 




r 

2Trr j&c 
4irr j 
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one dimensional 
cylindrical symmetry 
spherical symmetry 


(A-l) 


(A- 2) 


Here r x is a reference length large compared to a typical pore dimension, 
and A (r x ) is the total pore area normal to the mean flow at r = r , It 
is assumed that a is a constant property of the porous medium. This 
assumption requires that the rock be isotropic, and becomes more accurate 
with increasing 'values of r . 

The "average" velocity is defined 


v (r) 


E yr> 


JJ 


v (r,y,z) 

ir 


dA 

— p 


(A- 3) 


Here v is the "local" velocity in an individual pore. Similarly, the 
mass flow rate can be expressed 



1-26 


Q = [ I P (r,y,z) v (r,y,z) * dA 

J J P P P 

V rl 


= pv (r) A (r) (A-4) 

P 


where p is the local density in a pore. Now continuity equation can 
P 

be expressed, cf . Equation (1-5) , 


a p (r)v (r)r^ = Q (A-5) 


In order to satisfy Equation (A-5) , the average density is then defined as 


p (r) 


pv (r) 
v (r) 


(A-6) 


In the test of this report, the bars have been dropped from all 
average quantities. 



APPENDIX B* 
KNUDSEN NUMBER 


From the kinetic theory of gases , 


V = J b 2 XpV 


(B-l) 


where b 2 is a dimensionless constant of order one, and X is the mean 
free path. Invoking the perfect gas law, taking b 2 = 1, and 
rearranging Equation (B-l) , 


X = 2 UR T I 


V P 


(B-2) 


A Knudsen number based on L can be defined, Kn E 


_ X 


— . Then 
Xj 


Kn .asi 
L LV P 


(B-3) 


so that, for an isothermal flow, Kn , and p are inversely proportional, 

Li 


*As a general reference, see Chapman, S. and T. 0. Cowling (1964) , 

The Mathematical Theory of Non-Uniform Gases , Cambridge Univ. Press, 
New York, 



SYMBOLS 


a 


a i 


A = TTa 2 

A (r) 


A 

P 

A 

r 

b 


bi 


source radius 
constant (cf. Eq. 13) 
area of fluid source 
reference area 

total pore area normal to flow 

total cross-sectional area of rock 
constant [cf. Eq. (13) ] 

constant of order one appearing in slip boundary condition 


b 2 


Cl 


c 2 

C 

C 


a 

D 


D 


KA 


F 


j 


k 


constant of order one 

constant [cf . Eq. (13) ] 

constant [cf . Eq. (13) ] 

master diffusion coefficient 

dimensionless master diffusion coefficient 


1 

3 


X 

d 


average pore diameter, and tube diameter 

momentum diffusivity appropriate to a particular flow regime 
dV Knudsen diffusivity 

universal function of dimensionless pressure, equivalent to C 

index for one dimensional cylindrically symmetric, or 
spherically symmetry flow 

continuum permeability 

local Knudsen number based on pore (or tube) diameter 
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X 

i 

d 


±£_ 
tt yd 


/ 8RT 
/ TT 


reference Knudsen number based on pore (or tube) diameter 

local Knudsen number based on L 
tube length 

reference length for cylindrical ly symmetric flow 
square root of generalized permeability 
generalized permeability 
pressure 

background pressure 
reference pressure 
mass flow rate 

reference mass flow rate for symmetric flow 
length coordinate for symmetric flow 
reference length 

gas constant 

Reynolds number based on mass flow rate 
time 

temperature 
gas velocity 

local velocity in pore 

mean thermal speed of gas 
volumetric flow rate 
axial coordinate along tube 


area fraction 
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? = 


X 

X 
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- p = -X- dimensionless pressure 
2URT r Kn 

Xj 

mean free path 

mean free path at reference conditions 

viscosity 

density 

local density in pore 


T 

<P 

X 


4Dt 


a l* 


dimensionless time 
porosity 

weakly varying function of Kn^, assumed a constant of order one 



